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Abstract 

In the scenario of the electroweak baryogenesis we consider the dynamics of 
fermions with a spatially varying mass in presence of a CP-violating bubble wall 
and a uniform magnetic field perpendicular to the wall. The relevant quantity for 
baryogenesis, Rr-^l — Rl-*r , is studied {Rr^l and Rl^r being the reflection 
coefficients for right-handed and left-handed chiral fermions, respectively). 
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1. Introduction 



The connection between particle physics and cosmology has received much attention in 
the last years. In fact, cosmology can use the predictions of particle physics in order to 
solve cosmological problems. The origin of the asymmetry between matter-antimatter 
is one of those important open questions in cosmology. As Sakharov pointed out more 
than thirty years ago [T|, it is possible to start with a Universe in a baryon symmetry 
state and later dynamically to generate a net asymmetry by particle interactions. He 
postulated three conditions in order to explain the observed baryon excess: 1) Baryon 
number non- conservation; this violation occurs in the Standard Model through the 

axial anomaly, but the rate of baryon number non- conserving processes are exponentially 
suppressed at zero temperature. 2) C and CP-violation; this violation has been observed 
in the neutral K meson system and recently in the B meson decay jS] . 3) Departure from 
thermal equilibrium; in fact, in thermal equilibrium particles and antiparticles will have 
the same number density. 

The Standard Model has in principle all the necessary ingredients for Sakharov condi- 
tions for the generation of the baryon asymmetry of the Universe jl]. In particular the 
electroweak baryogenesis is a good candidate because give us an explanation in terms of 
experimentally accessible physics. Essentially the baryon asymmetry problem is twofold 
to understand: 1) why the observable density of baryons in the Universe is much greater 
than the density of the antibaryons and 2) why the density of baryons is much less than 
density of photons. Indeed, the measurement of the ratio of baryon density ub (defined 
as the number density of baryons minus the number density of antibaryons rij) and 
the photon density n^, at the present time, recently has been determined from the mea- 
surement of the relative heights of the first two CMB acoustic peaks by WMAP satellite. 
It takes the value 

5.9 X 10"^° < — < 6.4 X 10~^°. (1.1) 

CP-violation is crucial to generate the baryon asymmetry. It is well known, however, that 
the amount of CP-violation in the quark sector of the standard model cannot account for 
the observed baryon asymmetry [71IH]. In fact, the strength of CP- violation coming from 
the CKM phase is inadequate by 10-12 orders of magnitude {uB/n^ ~ 10~^°) [H]. 
Another question is connected with the third Sakharov condition. The baryogenesis nec- 
essarily requires non-equilibrium physics. In order to avoid the so called washout after the 
phase transition, it is necessary to have v{Tc)/Tc > 1, where v{Tc) is the vacuum expec- 
tation value of the broken Higgs field. This condition ensures the surviving of the baryon 
asymmetry after the walls has passed and indicates that the electroweak phase transition 
should be, if strong enough, of the first order |10j. Moreover, lattice studies have shown 
that for any Higgs mass, the phase transition would be so weak that sphaleron interactions 
are in equilibrium in the broken phase. This implies that the baryon asymmetry goes to 
zero immediately after its generation JT]. On the other hand, if we use one loop high 
temperature effective potential we have an Higgs boson mass Mh < 60 Gev ^2] which is 
in contrast with the experimental lower bound reached > 114, 4 Gev jT^]. Hence, the 
magnitude of CP-violation and the experimental limits on Higgs mass contradict the pos- 
sibility to have an asymmetry of the experimentally observed order within the framework 
of the Minimal Standard Model and one is lead to consider extensions of the Standard 
Model. If we demands that the electroweak phase transition has to be of first order, we 
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have a mechanism for the transition called bubble nucleation. It proceeds by formation 
and expansion of bubbles of new phase within the old ones. The study of the propagation 
of the bubble has been object of much attention as from the early works [14j, 15^. The 
first order provided a natural way to depart from equilibrium that takes places near the 
wall of the expanding bubbles. Moreover, in the unbroken phase where the sphaleron is 
active, the asymmetries of some local charges are converted into a baryon asymmetry. 
The baryon number flows into the new phase and the baryogenesis is produced. In this 
scenario of electroweak baryogenesis there are two interesting mechanisms that involve 
interacting fermions with CP-violating background of the bubble. In both cases, it is 
possible to have a net baryon flux by using fermion scattering off bubble walls ^Hj '■ The 
local baryogenesis jTTj, JH| in which the CP- violating and B- violating processes are near 
or at the bubble wall. 2) Non local baryogenesis PH] (or charge transport mechanism) in 
which CP-violation and baryon number violation are separated from one another. 
In our work we consider this last kind of electroweak baryogenesis. Moreover, the non 
local baryogenesis is effective for thin- wall regime jni-CH]- Hence, it is possible to treat 
as free the fermions that may only interact in a small region with Higgs background field. 
As a consequence of CP-violation on the wall, we will have asymmetric reflection and 
transmission, and the particles injected into the broken phase diffuse in the bubble. The 
correct way to treat the generation of the asymmetry in front of the bubble is controver- 
sial. A very simple mode is to consider an extra source of CP-violation, parameterized 
assigning a spatially varying complex mass to the Higgs field: 

m(t,x) = mR(t,x) + imj(t,x). (1.2) 

The imaginary part m/ gives rise to CP- violation. The problem to solve Dirac equation in 
the background of the bubble wall without CP-violation |'20'| and with CP-violation j2l] 
has been already studied. The authors also have studied the case of fermion scattering 
in the background of the bubble wall without CP-violation and in presence of a constant 
magnetic field perpendicular to the wall [22] ■ 

It is important to stress that the introduction of a magnetic field is not an academic 
exercise because astronomical observations show the presence of a cosmic magnetic field 
in all galaxies [22] • There are some possible mechanism for generating a cosmological 
magnetic field in the early Universe. One of these is related to the first order electroweak 
phase transition [21|, in which, bubbles of the new phase expanding in the old ones 
generate electric currents that, in turns, produce magnetic fields [2S|. 
In this paper we study and solve the Dirac equation in presence of a CP-violating planar 
bubble wall with a constant magnetic field perpendicular to the wall. We consider the 
scattering perpendicular to the wall and the radius of the bubble very large. Neglecting the 
time dependence of mass terms, we shall work in the hypothesis that m{t, x) only depends 
by the 2;-coordinate perpendicular to the wall: m{t, x) = m{z). The real part goes to zero 
when z —>■ —oo (symmetric phase) and goes to mo (broken phase) when z +oo, where 
mo is fermion mass. We calculate the reflection coefficients, Rr^l and Rl^r, of left- 
handed and right-handed fermions, respectively, at a radially expanding bubble. The 
difference AR = Rr^l ~ Rl^r is important as regards the baryon asymmetry [2ni in 
cosmology. 
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2. CP-violating Dirac equation 



In this Section we study tlie scattering of Dirac fermions off CP-violating bubble walls in 
the presence of an electromagnetic field A^. Thus, assuming that a fermion \E' is coupled 
to a complex scalar field $ through a Yukawa interaction with coupling gy, the Dirac 
equation reads: 



[ ip - m{t, x)Pr - m* {t, x)Pl - ] ^(t, f) = 0, 



(2.1) 



where e is the electric charge, Pr = (1 + 75)/2 and = (1 — 75)/2, are the right-handed 
and left-handed projection operators such that PrPl = PlPr = ^, Pr = Pr, Pi = Pl, 
and 

m{t,x) = -gYm,x)) (2.2) 

is a complex function of space-time. We only consider the motion of fermions perpendicu- 
lar to the wall, and the radius of the bubble very large in order to approximate the surface 
by a fiat wall perpendicular to the z-axis. Hence, we shall consider m{t,x) = m{z). In 
order to solve Eq. ()2.H) . we make the following ansatz: 



(2.3) 



where a = ±1 for positive and negative energy solutions, respectively. Moreover, we 
assume that corresponds to a constant and uniform magnetic field directed along the 
2;-axis with strength B. We can then choose the Landau gauge 



= (0,0,-Px,0). 
Inserting Eq. (Q into Eq. (EU, we get 



E^ + — -\m\^ + ij 

dz"^ dz 



ndrriR o dmi . „ i 9 



dz 



0. 



(2.4) 



(2.5) 



It is useful to introduce a parameter mass term a, (1/a being the characteristic size of the 
thickness of the wall) and define 



mR{z) = mof{az) = mof{x), 
mi{z) = mog{az) = mog{x), 



(2.6) 
(2.7) 



where 



X = az, T = at, e = E/a, ^ = mo/a, h = eB/a?. (2-8) 
Taking into account Eqs. ()2.6|) - ()2.8p . the ansatz ()2.3|) and Eq. ()2.5|) become, respectively 



^lir^x) 



ae-f° + «7^-r- + ^/ - ^^5-75 + bx-/'^ 



dx 



e{f + 9')+iil'r-i 757'^' + 



3^' 



,1.2 



dx"^ 



0, 



(2.9) 
(2.10) 



where a prime denotes differentiation with respect to x. Let us suppose that \g{x)\ ^ 1. 
The small value of g allows to consider it as a perturbation and to consider only first 
order terms in g. We assume that 



lim f{x) = 1, lim f{x) = 0. 



x— >+oo 



x—>—oo 



(2.11) 
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In first order approximation we write 

The wave function ip^^^ is a solution of the unperturbed equation 

V'^^^x) = 0, 



(2.12) 



(2.13) 



which is obtained from Eq. (j2.1()|l with g = 0. The perturbation ip^^^ can be calculated as 



^(1) (x) = I dx' G{x, x')V{x') (x') 



(2.14) 

where V{x) = — ^fi'(x')757^, and G{x,x') is a Green function that satisfies the following 
equation, 



dx'^ 



Gi3^{x,x') = —6a-yS{x — x'). (2.15) 



J a/3 



Inserting Eq. (jTT^ in Eq. (jTTIH) we have 



^('r,x) 



d 



■,(0) 



(2.16) 



to the first order in g. In order to obtain ^^^^ we have to solve Eq. ()2.13|) (see Refs. 
). Let us expand ip^'^^ in terms of the eigenstates of 7'^: 



^(o)(a;) =<p]^>{x)u 



4, 



(2.17) 



where s = 1, 2, and 



1 
72 



/ 1 \ 



±i 
\ / 



V2 



( \ 

1 




f2.18) 



The spinors satisfy the relations 



7°«± = n% , 
7V< = 



7^M± = ±m 



± 1 



75M± = T«(-1)X' 
that will be useful in the following. Inserting Eq. ()2.17p into Eq. ()2.13|) . we get 



dx^ 



4'^(x) = 0. 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



(2.23) 
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We indicate the independent solutions of Eq. ()2.23|) as (p^±"'\x) and ^^±'^'\x). The 



asymptotic properties of (f)^^"'''^ and are, respectively 



(x) 



for X +00, 
for a; —>■ —00, 



-OsX 



for X +00, 
g-/3s2: for X —00, 



where 



(3s = i^e^-{-iyh, 



(2.24) 
(2.25) 



(2.26) 
(2.27) 



and 7±(as,/3s) are constants such that 7±(q!s,/9s)* = 7±(— Os, — /9s). The general unper- 
turbed solution is 



(2.28) 



=1,2 



The incident wave function coming from x = — cxo is reflected and transmitted, while at 
X = +00 there is only the transmitted wave. Thus, we have A^""^) = for a = +1 and 
= for a = -1. 

Following Ref. [21], in order to calculate the Green function G'(x,x'), we introduce a 
unitary matrix 

/II \ 



U ={u\ u\uiul) = ^ 



V2 



11 

1 
\ -i i / 



(2.29) 



and we write G(x, x') as 

/ G_i_(x,x') 

G(x, x') = U 



GJx.x') 



G+(x, x'l 



v 



G_(x, x') J 



(2.30) 



where 



dx"^ 



G±{x, x') = — 5(x — x'). 



(2.31) 



Following the standard method [27j, we find for the Green function the expression 



G^'"Vx,x') 



__£_[0(:-."^)(a;) + 4^''^)0(+'^'"^)(x)] 0i-"'"^)(x') if X < 



X 



(2.32) 



X' + C, 



(X')]0lr"'"^^(x) ifx>x'. 
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Therefore, if we know the Green function and wave function ■ip^'^^ it is possible to find, by 
Eq. (j2.14j) . the wave function ip^^^: 



(+) /N A+as), 



+ 



x) + c 



(+UL+"^)(a;) / dx'g'{x')<j)'-^'''\x')<P^^'''\x')}. (2.33) 



Let us introduce the following quantities: 



+00 



r(^) 



00 
+00 



dxg\x) 

dxg\x) (t)t"^\x). 



(2.34) 
(2.35) 



In Appendix A we shall calculate the asymptotic expression of the transmitted, incident 
and reflected wave functions (for definiteness, let us consider the case a = +1). They are 
respectively. 



X U^-tt 



leT+OsX 



+ e 



(-ime+a^ 



2eas 



2eas{( - as) 



2e 



2ea, 



2e 



u 



(2.36) 



[^siT,x)fZ_,, = 4+)7+(«.,/?.)e-^^^+^^^ 



X <-/3s 



-l)^ee^N-(-«s,/5.) , ^bc^^h^'^ 



+ e 



2asPs 7+(«s 
r(^) 



+ 



2eas{^ — a 



2 ^ i-mg- 



2e 



1 + 



lymr 7-(-«s,/5s) i-mg- 



2ea. 



2e 



u 



[^s{T,x)l%, = [^s{r,x)y:: 



+1 



(2.37) 
(2.38) 



where g± = \imx^±oog{x). From Eqs. ()2.3(i|) - ()2.38|) we calculate the vectorial current 
jy = il'7^\l/ and the axial current = '^'j'^'j^'^. After some manipulations we obtain 



ifyj^^^ = 2e|AW|>,|(l + 5*-'^), 
{fyX^ = 2e|AWn7+(a„/3,)n/3,|(l + 5r + 5r) 



f -3 \ refl /' „■ 3 ^ inc 

{jV,s) - 



iJvJ 



(2.39) 
(2.40) 
(2.41) 
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where 



^tran 



2e\a\ 
(-1)' 



-Re 



e a. 



mc 
2 



2e|«,||A 

The transmission and reflection coefficients are 



+ Im 



,(+) 



r(^) 



Ts,B 
Rs,B 



( a3 \tran 

Uv,s) 

_(je>r 



T, 



(0) 



s,B 



(0) 
s,B 



1 I rrcfl I rrefi rinc rmcN 

1 + di +62 -di -62 ), 



where 



T. 



(0) 



s,B 





as 




m 


l7+(" 




2 



R 



(0) 



l+ias,Ps) 



(2.42) 
(2.43) 

(2.44) 

(2.45) 
(2.46) 

(2.47) 



are the transmission and reflection coefficients in the absence of CP- violation [221, 

5f = -4-1 (2.48) 



rrcfl ein 
^1 = 



The unitary condition, + T^'^ = 1, holds as it should be. We find the expression 
for ci^^ by requiring the unitary condition for the refiection and transmission coefficients 
Rs^B and Ts^b- Imposing that 

Rs,B + Ts,B = l, (2.49) 

we get 

M^L. (2.50) 
i Re|4''/(« - a.)] 

Now, we calculate the axial currents (j^^)™'^ and {jAsy^^^ analogously to the vectorial 
current. We get the following expressions: 

X [e^{l + 6r) + \(3sWl + 6r + 26'^'')], (2.51) 



/ -3 Nrcfl _ / -3 Nine | 

where, here and in the following (see Eqs. ()2.58j) . ()2.(i3|l and ()2.fi4j) ). we set g- - 
Finally, we are able to calculate the relevant quantity for the generation of 
logical baryon asymmetry, 



where 



R 



\refl 



,refl 



,refl 



R 



L^R 



( riS \ 'mc 

Ur,s) 



Ul,s) Ua,s) Uv,s) 



(2.52) 

= 0. 

a cosmo- 
(2.53) 

(2.54) 
(2.55) 



Taking into account the expressions for the currents, we obtain 
where we have defined the quantities 



rmc 



2\as\ [ -f+{as,(3s) 



c.c. 



rrefi rinc 



(2.56) 



(2.57) 



In the case of vanishing magnetic field, AR^'^^ reduces to the same expression calculated 
in Ref. EH- 

In order to study AR^^^ we have to know f(x) and g{x). As regards f{x) we take the 
usual bubble wall profile 

„ , , 1 + tanhx 

fix) = ^ . (2.58) 

With the profile ()2.58|) . the solutions (j)^^"''\x) of the unperturbed equation ()2.23|) may 
be written as 



y 2 (1 - y) 2 2F1 



7; ^ 1> 7^ , l-as]y 



(2.59) 



{y) = y"'(i -y) 2F1 



7; \- 1, 7; , l + as;y 



(2.60) 



where y = 1 — /(x), and 2Fi{a, b, c; x) is the the hypergeometric function. Now, 7±(q;s, Ps 
and R^B explicitly given by 

and 



r(-a, + 1) Ti-P, 



R 



T[{-as -f3s± 0/2] r[(-«, -(3.,T 0/2 + 1] 
(0) _ sin[(7r/2)(a, - + 0] sin[(7r/2)(a, - - 0] 



(2.61) 



(2.62) 



^'^ sin[(7r/2)(«, + Ps + 0] sin[(7r/2)(a, + - 0] ' 
The functional form of g{x) is unknown. Following Ref. [21], we consider two expres- 



sions for g{x): 



g{x) = Aefix) 
g{x) = AOfix) 



(2.63) 
(2.64) 



where the parameter characterizes the magnitude of CP-violation. 

In Figs. 1-4 we plot AR^^'' / A9 as a function of normalized energy e* = ea* = E/rrio, 
at fixed thickness a* = a/rrio = 1/^ [a* = 5 in Figs. 1 and 3, a* = 1 in Fig. 2 and 4) for 
different values of the magnetic field b* = ba*"^ = eB /m^. We take s = 1 (left panels) and 
s = 2 (right panels), with g{x) = A9p{x) in Figs. 1,2, and g{x) = A6f\x) in Figs. 3,4. 
We found that AR^^^ display the following peculiar properties: 

i) At fixed thickness of the wall (proportional to 1/a*) and magnetic field, the absolute 
value of Ai?(^) goes to zero when the energy of the incident particles approaches to infinity; 

m) The maximum value of |Ai?(^)| vary both with the thickness of the wall and with 
the functional form of g{x). In particular, we found that at fixed energy and magnetic 
field, |Ai?^*^| is an increasing function of the thickness of the wall; 
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in) The presence of a magnetic field generates a reflection asymmetry between spin-up 
and spin-down particles. In particular, the effect of the magnetic field is to shift the values 
of Ai?(*) (with respect to the case B = 0) towards lower energies in the case s = 1, and 
higher energies in the case s — 2. 

The peculiar global quantity in non-local defect mediated electroweak baryogenesis is 
the flux of Icpton number radiated by the bubble wall. In the thermal frame of reference, 
it is given by (see e.g. Ref. [16]): 




(2.65) 



where T is the temperature and, considering only one type of particle, L{k) = /|i?(k)p. 
Here, R(k) is the reflection amplitude for a particle of momentum k, is the angle of 
reflection off the advancing wall, k = |k|, and I is the lepton number. In the case of an 
infinitely planar wall only the motion of fermions perpendicular to the wall is important, 
and then we have $l(-B) oc J dE e~^^^AR{E, B) , where E is the energy of the scattered 
particle, and Ai? = Rr^l — Rl^r ■ Indeed, the relevant quantity for the generation of a 
cosmological baryon asymmetry turns out to be $^°*(5) -|- $^*(— 5), where 

$f (S) = ^ ^t\B) ^otY, I dEe-^/^ AR('\E,B) . (2.66) 

s=l,2 8=1,2'^ 

Here, a > is a numerical factor, and s refers to the spin of the scattered particle. It is 
straightforward to verify that the following transformation law hold: 

AR^'\-B) = AR^'\B), (2.67) 

where s = 1 if s = 2, and s = 2 if s = 1. Finally, taking into account the previous 
equations, we have 

^i\B) + ^i\-B) = 4^)(S) + $f (S) = $f (S), (2.68) 

and 

¥£\B) + $f (-S) = 2$f (S). (2.69) 

In general, the quantity ^*£^[B) is not null. For example, from Fig. 1, we get that ^^1\b) 
and ^^1\b) are strictly negative quantities. Hence, in this case we have $^°*(i?) = 
^^1\b) + ^^1\b) < 0, and then, in general, the total lepton number flux radiated by an 
inflnitely planar wall is different than zero. 

3. Conclusions 

In this work we have investigate a baryogenesis scheme at a flrst order electroweak phase 
transition in which fermions interact with a CP- violating thick bubble wall in the back- 
ground of an uniform magnetic field perpendicular to the wall. 

We have studied the net fiux by the difference AR^'^^ of the refiection coefficients of the 
left-handed and right-handed fermions incident from the unbroken phase at the expanding 
bubble wall. We have found that the presence of a magnetic field generates a reflection 
asymmetry between spin-up and spin-down particles. Moreover, we have showed that the 
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only effect of the magnetic field is to shift the values of AR^^^ (with respect to the case 
of vanishing magnetic field) towards lower energies in the case of spin-up particles, and 
higher energies in the case of spin-down particles. 

In non-local defect mediated electroweak baryogenesis, the total baryon asymmetry 
turns out to be proportional to the total lepton number flux radiated by a bubble wall, 
— $^*(— S). We have showed that, in general, this quantity is not null. In any 
case, the cosmological implications of a non-vanishing lepton number flux are beyond the 
aim of the present paper. We deserve such an analysis to a future work. 
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Appendix 



A Transmitted wave function 



In this Appendix we calculate the transmitted wave function. In the same way it is 

straightforward to calculate the incident and reflected wave functions. 

The wave function to the first order is given by Eq. ()2.1fi|) . that we rewrite for simplicity: 



(Al) 



Hence, in order to calculate the asymptotic form of \E'(r, x), it is necessary to know the 
asymptotic form of ip^^\x) and ^^-^\x). Taking into account Eqs. ()2.24|) and ()2.28|) we 
obtain, for positive energy solutions and in the limit x +00, 



(A2) 



As regards ^^^\x) we start with Eq. ()2.33|) . Taking into account Eqs. ()2.24|) . ()2.25|) . and 
Eqs. (ESI, (1^^ we have 



-m 



2a. 



u e 



(A.3) 



Remembering that lim^^^^ f {x) = 1, the transmitted wave function turns out to be 



(A.4) 



Finally, inserting Eqs. ()A.2|) and ()A.3|) into Eq. ()A.4|) . we get 



+ 



X e 



2a 

-itT+asX 



(A5) 



Following the same arguments of Appendix in Ref. [2T], it is possible to redefine A^^"^ as 



-1)'— ^- 

^ 2e^ 



(A.6) 



Taking into account Eq. ()A.6|1 . it is straightforward to cast Eq. ()A.5jl into the form of 

Eq. (ESSD- 
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Figure 1: Ai?*^*Y^^ versus e*, with s = 1 (left panel) and s = 2 (right panel), for different 
values of normalized magnetic field h* , in the case g = ^dp and a* = 5. Thick solid line: 
b* — 3; thick long dashed line: b* = 2; thick short dashed line: b* = 1.125; thin solid line: 
b* — 1; thin long dashed line: b* —0.5; thin short dashed line: fo* = 0. 
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Figure 2: AR^^'' / A9 versus e*, with s = 1 (left panel) and s = 2 (right panel), for different 
values of normalized magnetic field b*, in the case g = A9p and a* = 1. Thick solid line: 
b* = 1.2; thick long dashed line: b* = 1.1; thick short dashed line: b* = 1.02; thin solid 
line: b* = 1; thin long dashed line: b* = 0.5; thin short dashed line: b* = 0. 
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Figure 3: AR^'^/A9 versus e*, with s = 1 (left panel) and s = 2 (right panel), in the case 
g — A9f' and a* — 5. The values of b* are the same as in Fig. 1. 
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Figure 4: AR^'^Ae versus e*, with s = 1 (left panel) and s — 2 (right panel), in the case 
g — AOf and a* — 1. The values of h* are the same as in Fig. 2. 
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